arXiv:1504.06436vl [cond-mat.mes-hall] 24 Apr 2015 


Second harmonic coherent driving of a spin qubit in a Si/SiGe quantum dot 

P. Scarlino^, E. Kawakami^, D. R. Ward^, D. E. Savage^, M. G. Lagally^, Mark 
Eriesen^, S. N. Coppersmith^, M. A. Eriksson^ and L. M. K. Vandersypen^ 

^ Kavli Institute of Nanoscience, TU Delft, 

Lorentzweg 1, 2628 CJ Delft, The Netherlands 
^ University of Wisconsin-Madison, 

Madison, WI 53706, USA 

(Dated: April 27, 2015) 

We demonstrate coherent driving of a single electron spin using second harmonic excitation in a 
Si/SiGe quantum dot. Our estimates suggest that the anharmonic dot confining potential combined 
with a gradient in the transverse magnetic field dominates the second harmonic response. As 
expected, the Rabi frequency depends quadratically on the driving amplitude and the periodicity 
with respect to the phase of the drive is twice that of the fundamental harmonic. The maximum 
Rabi frequency observed for the second harmonic is just a factor of two lower than that achieved 
for the first harmonic when driving at the same power. Combined with the lower demands on 
microwave circuitry when operating at half the qubit frequency, these observations indicate that 
second harmonic driving can be a useful technique for future quantum computation architectures. 

PACS numbers: 73.21.La, 71.70.Ej, 72.25.Rb, 75.70.Tj 


Gontrolled two-level quantum systems are essential el¬ 
ements for quantum information processing. A natural 
and archetypical controlled two-level system is the elec¬ 
tron spin doublet in the presence of an external static 
magnetic field [Ei- The common method for driving 
transitions between the two spin states is magnetic res¬ 
onance, whereby an a.c. magnetic field (Ba.c.) is applied 
transverse to the static magnetic field with a fre¬ 

quency, fMWi matching the spin Larmor precession fre¬ 
quency /l = gUBBtotjh {h is Planck’s constant, is the 
Bohr magneton and Btot the total magnetic field acting 
on the spin). When the driving rate is sufficiently strong 
compared to the dephasing rates, coherent Rabi oscilla¬ 
tions between the ground and excited state are observed. 


Both spin transitions and Rabi oscillations can be 
driven not just at the fundamental harmonic but also at 
higher harmonics; i.e., where the frequency of the trans¬ 
verse a.c. field is an integer fraction of the Larmor fre¬ 
quency, fMW = Il/b, with n an integer. Second or 
higher harmonic generation involves non-linear phenom¬ 
ena. Such processes are well known and explored in quan¬ 
tum optics using non-linear crystals 0 and their selec¬ 
tivity for specific transitions is exploited in spectroscopy 
and microscopy 0^- Two-photon transitions have been 
extensively explored also in superconducting qubit sys¬ 
tems [i-liil- In cavity QED systems, a two-photon pro¬ 
cess has the advantage that it allows the direct transition 
from the ground state to the second excited state, which 
is forbidden in the dipole transition by the selection rules 


Eor electron spin qubits, it has been predicted that 
the non-linear dependence of the ^-tensor on applied 
electric fields should allow electric-dipole spin resonance 


(EDSR) at subharmonics of the Larmor frequency for 
hydrogenic donors in a semiconductor [l3, [lH . Eor elec¬ 
trically driven spin qubits confined in a (double) quan¬ 
tum dot, higher-harmonic driving has been proposed that 
takes advantage of an anharmonic dot confining potential 
[iMi or a spatially inhomogeneous magnetic field (2^ . 
Although higher harmonics have been used in continu¬ 
ous wave (GW) spectroscopy for quantum dots hosted 
in GaAs, In As, InSb and carbon nanotubes 2ll-[2^. co¬ 
herent spin manipulation using higher harmonics has not 
been demonstrated previously. 

In this letter we present experimental evidence of co¬ 
herent second harmonic manipulation of an electron spin 
confined in a single quantum dot (QD) hosted in Si/SiGe 
quantum well. We show that this second-harmonic driv¬ 
ing can be used for universal spin control, and we use 
it to measure the free-induction and Hahn-echo decay of 
the electron spin. Eurthermore, we study how the sec¬ 
ond harmonic response varies with the microwave am¬ 
plitude and phase, and comment on the nature of the 
non-linearity that mediates the second harmonic driving 
process in this system. 

The quantum dot is electrostatically induced in an un¬ 
doped Si/SiGe quantum well structure, through a combi¬ 
nation of accumulation and depletion gates (see Sec. I of 
[m for full details). The sample and the settings are the 
same as those used in Ref. [27|. A cobalt micromagnet 
next to the QD creates a local magnetic field gradient, 
enabling spin transitions to be driven by electric fields 


All measurements shown here are performed using 
single-shot readout via a QD charge sensor [s^. They 
make use of four-stage gate voltage pulses implement- 
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FIG. 1. (a) Measured resonance frequencies as a function of 
externally applied magnetic field Bext- The long microwave 
burst time tp =700 /iS ^ T 2 means that the applied excitation 
is effectively continuous wave (CW). The microwave source 
output power was P — —33 dBm to —10 dBm (—20 dBm to 
—5 dBm) for the case of fundamental (second) harmonic ex¬ 
citation, decreasing for lower microwave frequency in order to 
avoid power broadening. The red and green lines represent fits 
with the relation hf = {Bext — B\\)^ + respectively 

to the resonance data labeled (2) and (3) (we excluded points 
with Bext < 700 mT from the fit because the micromagnet 
apparently begins to demagnetize there) [^. (b) Schematic 
of the energy levels involved in the excitation process, as a 
function of the total magnetic field at the electron location. 
The dashed arrows correspond to the four transitions in panel 
(a), using the same color code, (c) Schematic of an anhar- 
monic confinement potential, leading to higher harmonics in 
the electron oscillatory motion in response to a sinusoidally 
varying excitation, (d) Measured spin-up probability, P^, as 
a function of applied microwave frequency, /mpv, for Bext — 
560.783 mT (P = —30 dBm for the fundamental response, 
P = —12 dBm for the second harmonics), averaged over 150 
repetitions per point times 80 repeated frequency sweeps (160 
mins in total). The frequency axis (in red on top) has been 
stretched by a factor of two for the second harmonic spin re¬ 
sponse (red datapoints). From the linewidths, we extract a 
lower bound for the dephasing time — 760 ± 100 ns, 

= 810 ± 50 ns, = 750 ± 40 ns and = 910 ± 
80 ns. The Gaussian fits through the four peaks use the same 
color code as in panels (a) and (b). 


ing (1) initialization to spin-down, (2) spin manipu¬ 
lation through all-electrical microwave excitation, (3) 
single-shot spin readout, and (4) a compensation/empty 
stage [13 . The results of many single-shot cycles are used 
to determine the spin-up probability, P|^, at the end of 
the manipulation stage. 


First we apply long, low-power microwave excitation 
to perform quasi-CW spectroscopy. Fig. [TJa) shows four 
observed spin resonance frequencies, through as 
a function of the external magnetic field. The resonances 
labeled (1) and (2) represent the response at the fun¬ 
damental frequency. As in [271, these two closely spaced 
resonances correspond to the electron occupying either of 


the two lowest valley states, both of which are thermally 
populated here. The other two resonances occur at ex¬ 
actly half the frequency of the first two, = 2/q^\ 
= 2/q^\ and represent the second harmonic re¬ 
sponse. 

The effective ^f-factors extracted from the slopes for the 
second harmonic response are half those for the first har¬ 
monic response [see Fig. [TJa) inset]. The relevant energy 
levels and transitions as a function of the total magnetic 
field, Btoti are visualized in Fig. [IJb), where the color 
scheme used for the resonances is the same as in Fig.[Tl[a). 
We see two sets of Zeeman split doublets, separated by 
the splitting between the two lowest-energy valleys, Ey. 
The transition between the Zeeman sublevels within each 
doublet can be driven by absorbing a single photon or two 
photons, as indicated by the single and double arrows. 

To drive a transition using the second harmonic, a non¬ 
linearity is required. In principle, several mechanisms 
can introduce such a non-linearity in this system (see 
Sec. II of [HI). First, as schematically shown in Fig.[T][c), 
if the confining potential is anharmonic, an oscillating 
electric field of amplitude Ea.c. and angular frequency 
uj = 271 fMW induces effective displacements of the elec¬ 
tron wavefunction with spectral components at angular 
frequencies ncj, with n an integer. In analogy with non¬ 
linear optical elements, we can look at this process as 
generated by an effective non-linear susceptibility of the 
electron bounded to the anharmonic QD confinement po¬ 
tential. 

The gradient in the transverse magnetic field in the 
dot region {B± in green) converts the electron motion 
into an oscillating transverse magnetic field of the form 

5^-^-(t) = B^ cos(cjt) + B 2 uj cos(2cjt) + ... (1) 


that can drive the electron spin for huj = Ez^ 2huj = Ez 
and so forth 0- A possible source of anharmonicity 
arises from the nonlinear dependence of the dipole mo¬ 
ment between the valle y ( or valley-orbit) ground (v-) 
and excited states (i;+) [3l|, as a function of F^a.c.- 

A second possible source of nonlinearity is a variation 
of the transverse field gradient, with position [see 

Fig. me)]. Even if the confining potential were harmonic, 
this would still lead to an effective transverse magnetic 
field containing higher harmonics, of the same form as 
EqlU 

A third possibility is that not only the transverse mag¬ 
netic field but also the longitudinal magnetic field varies 
with position. This leads to an a.c. magnetic field which 
is not strictly perpendicular to the static field, which is in 
itself sufficient to allow second harmonic driving 32, 33| . 
even when the confining potential is harmonic and the 
field gradients are constant over the entire range of the 
electron motion. 

However, simple estimates indicate that the second 
and third mechanisms are not sufficiently strong in the 
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FIG. 2. Rabi oscillations, (a) Measured spin-up probability, 
Pf, as a function of microwave burst time {Bext = 560.783 
mT, fMW = 6.4455 GHz) at four different microwave pow¬ 
ers, corresponding to a rms voltage at the source of 998.8 
mV, 1257.4 mV, 1410.9 mV, 1583.0 mV. (b) Rabi frequencies 
recorded at the fundamental harmonic, (blue triangles, 
adapted from [13]), and at the second harmonic, (green 
squares), as a function of the microwave amplitude emitted 
from the source (top axis shows the corresponding power). 
For the second harmonic, the amplitude shown corresponds 
to a 5 dB higher power than the actual output power, to com¬ 
pensate for the 5 dB lower attenuation of the transmission line 
at 6 GHz versus 12 GHz (estimated by measuring the coax 
transmission at room temperature). The green solid (dashed 
black) line is a fit of the second harmonic data with the rela¬ 
tion log(/^) oc 21og(F;a.c.) [log(/^) oc log(F;a.c.)]. The large 
error bars in the FFT of the data in Fig.[2lja) arise because we 
perform the FFT on only a few oscillations. Bext— 560.783 
mT. 
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FIG. 3. Phase control of oscillations, (a) Probability mea¬ 
sured after applying two 7r/2 rotations via second harmonic 
excitation, as a function of the relative phase between the 
two microwave bursts, Acj). The two rotations are separated 
by T = 100 ns (black) and t — 2 jis (red). (P = 16.0 dBm, 
= 560.783 mT, fnw = = 6.44289 GHz), (b) Sim- 

ilar to panel (a), but now driving the fundamental harmonic 
for T = 20 ns (black) and t — 2 jis (red). (P = 12.0 dBm, 
Bext = 560.783 mT, fuw = = 12.88577 GHz). Inset: 

Microwave pulse scheme used for this measurement, (c) Mea¬ 
sured spin-up probability, P^ (1000 repetitions for each point), 
as a function of /mvc and the relative phase A0 between two 
7r/2 microwave bursts (130 ns, P= 16.0 dBm) for second har¬ 
monic excitation, with r = 50 ns. The measurement extends 
over more than 15 hours. 


present sample to allow the coherent spin manipulation 
we report below (see Sec. II of [HI). We propose that the 
first mechanism is dominant in this sample, supported by 
our observation that the strength of the second harmonic 
response is sensitive to the gate voltages defining the dot. 


In Fig. Hd) we zoom in on the four CW spin res¬ 
onance peaks, recorded at low enough power to avoid 
power broadening (see Sec. I of [Hj). Fitting those res¬ 
onances with Gaussians, we extract the dephasing times 
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-, giving val- 


ues in the range of 750 to 910ns for all four resonances 
[see caption of Fig. [Il[d)]. This directly shows that the 
linewidth (FWHM) extracted for the two-photon pro¬ 
cess is half that for the one-photon process, as expected 

11011 . 


From the relative peak heights in Fig. Ed) , we can 
estimate the ratio of the Rabi frequencies between the 
two peaks in each pair (see Sec. I of [HI). In we 
found that the relative thermal populations of the two 
valleys (e^^^/e^^^) were about 0.3 ±0.1 to 0.7±0.1. Given 
this, the ratio between the Rabi frequencies, /i, extracted 
from the peak heights is rR{2ph) = i 0-2 

for the second harmonics. This is different from the 
ratio observed in for the fundamental harmonic, 
rnilph) = = 1-70 ±0.05 0. 

Such a difference is to be expected. The ratio rii{2ph) 


is affected by how the degree of anharmonicity in the 
confining potential differs between the two valleys. In 
contrast, rR{lph) depends on how the electrical suscep¬ 
tibility differs between the two valleys [35[. In addition, 
since the valleys have different charge distributions (H| . 
the microwave electric field couples differently to the two 
valley states, and this difference can be frequency depen¬ 
dent 00. Because the second harmonic Rabi oscilla¬ 
tions are driven at half the frequency of the Rabi oscil¬ 
lations driven at the fundamental, this frequency depen¬ 
dence also contributes to a difference between rR{lph) 
and rR(2ph). We note that the difference in Rabi fre¬ 
quency ratio between the 1-photon and 2-photon case 
demonstrates that the second harmonic response is not 
just the result of a classical up-conversion of the mi¬ 
crowave frequency taking place before the microwave ra¬ 
diation impinges on the dot, but takes place at the dot 
itself. 

The second harmonic response also permits coherent 
driving, for which a characteristic power dependence is 
expected HI, S, [HI . Fig. [2][a) shows Rabi oscillations, 
where the microwave burst time is varied keeping /mw = 
/q ^ for different microwave powers. We note that the 
contribution to the measured spin-up oscillations coming 
from the other resonance, (4), is negligible because the 
respective spin Larmor frequencies are off-resonance by 
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2 MHz, ~ 1 population is ^ three times 

smaller. 

To analyze the dependence of the Rabi frequency on 
microwave power, we perform a FFT of various sets of 
Rabi oscillations similar to those in Fig. [2{a). Fig. [Sfb) 
shows the Rabi frequency thus obtained versus mi¬ 
crowave power for driving both at the second harmonic 
(green) and at the fundamental (blue), taken for identical 
dot settings . We see that for driving at the frequency 
of the second harmonic, the Rabi frequency is quadratic 
in the applied electric field amplitude (linear in power), 
as expected from theory [20|, l33|, l38| . When driving at 
the fundamental resonance, the Rabi frequency is linear 
in the driving amplitude, as usual. It is worth noting 
that at the highest power used in this experiment the 
Rabi frequency obtained from driving the fundamental 
valley-orbit ground state spin resonance is just a factor 
of two higher than the one from driving at the second 
harmonic. This ratio indicates that the use of second 
harmonic driving is quite efficient in our device. 

A further peculiarity in coherent driving using second 
harmonics is seen when we vary the phase of two con¬ 
secutive microwave bursts. Fig. [H^a) shows the spin-up 
probability following two 7r/2 microwave bursts with rel- 
ative phase A0, resonant with /q ’ and separated by 
a fixed waiting time r. For short r, the signal oscil¬ 
lates sinusoidally in A0 with a period that is half that 
for the single-photon case [compare the black traces in 
Figs.[Ha,b)]. 

Therefore, in order to rotate the electron spin around 
an axis in the rotating frame rotated by 90 degrees with 
respect to a prior spin rotation axis (e.g. a Y rotation 
following an X rotation in the rotating frame), we need to 
set A0 to 45 degrees, instead of 90 degrees, when driving 
via the second harmonic. Of course, for r ^ T 2 *, the con¬ 
trast has vanished, indicating that all phase information 
is lost during the waiting time [Fig. [3][a,b) red traces]. 
Fig. El^c) shows two-pulse measurements as in Fig. Ha) 
as a function of frequency detuning and phase difference, 
where we can appreciate the extraordinary stability of 
the undoped device. 

To probe further the coherence properties of the spin, 
we perform a free induction (Ramsey) decay measure¬ 
ment, see Fig. ^ function of frequency detuning 

and delay time, r, between the two bursts. The absence 
of a central frequency symmetry axis is due to the pres¬ 
ence of two superimposed oscillating patterns, originat¬ 
ing from the resonances at and Figs. H^-e) 

show sections of the Ramsey measurement in Fig. 
corresponding to different waiting times r (see the white 
dashed lines). The visibility of the Ramsey fringes clearly 
decreases for longer waiting times between the two 7r/2 
pulses. Fitting the decay of the visibility of the fringes 
as a function of r with a Gaussian (oc exp[—(t/T 2 *)^], see 
Sec. I of [HI) we find T 2 = 780 ± 110 ns, in agreement 
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FIG. 4. Ramsey fringes, (a) Measured spin-up probability, 
as a function of Jmw and waiting time r {Bext = 560.783 
mT, P = 13.0 dBm) between two 7r/2 pulses (130 ns) with 
equal phase, showing Ramsey interference. Each data point is 
averaged over 300 cycles. Inset: Microwave pulse scheme used 
for this measurement, (b) Fourier transform over the waiting 
time, T, of the data in panel (a), showing a linear dependence 
on the microwave frequency, with vertex at /mw = and 
slope fRamsey = 2A/Mw(black dashed lines). The expected 
position of the FFT of the signal arising from resonance is 
indicated by the dotted black line. For comparison, the white 
dashed line represents the relation fRamsey = A/mw- (c-d- 
e) Sections of the Ramsey interference pattern in (a) along 
the three white dashed lines; the respective waiting times are 
indicated also in the inset of each panel, (f) Measured spin-up 
probability as a function of the total free evolution time, r, in 
a Hahn echo experiment (pulse scheme shown in inset). The 
decay curve is fit well to a single exponential (blue). Here, 
fuw = Be^t= 560.783 mT. 


with the value extracted from the linewidth. 

Furthermore, and analogously to the observations of 
Fig. H^), we report a doubling effect in the frequency 
of the Ramsey oscillations, fRamseys as a function of 

the detuning A/mvv = Imw - ■ Fig. m shows 

f Ramsey extracted from the data in Fig. IDJa) 

via a FFT over the waiting time r. The black dashed 
line indicates the condition fRamsey = closely 

overlapping with the position of the yellow peaks in 
the FFT. The black dotted line indicates the condition 
f Ramsey = “ /o^4 f^is second resonance is not 

very visible in the data, due to the lower population of the 
corresponding valley. For comparison, the white dashed 
line indicates the condition fRamsey = which is 

the expected response when driving at the fundamental. 
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Finally, we perform a spin echo experiment via second 
harmonic driving. Fig. SJ^f) shows as a function of 
the total free evolution time r, for a typical Hahn echo 
pulse sequence (illustrated in the inset) consisting of 7r/2, 
TT and 7r/2 pulses applied along the same axis, separated 
by waiting times r /2 j39|. A fit to a single exponential 
yields = 20.6 ± 6.5/is at Bext= 560.783 mT, com¬ 

patible with the of 23.0 ±1.2 /rs we observed at the 
same magnetic field when driving via the fundamental 
harmonic (see Sec. I of [HI). 

To summarize, we report coherent second harmonic 
driving of an electron spin qubit defined in a Si/SiGe 
quantum dot, including universal single-spin rotations. 
The non-linearity that permits second harmonic driving 
is likely related to the anharmonic confining potential for 
the electron. This means that routine use of second har¬ 
monics for spin control would be possible provided there 
is sufficient control over the degree of anharmonicity of 
the confining potential. This could be very useful since 
driving a spin qubit at half its Larmor frequency would 
substantially simplify the microwave engineering required 
for high fidelity qubit control. 

We acknowledge M. Rudner, A. Palyi and T. Jullien 
for useful discussions, R. Schouten and M. J. Tiggelman 
for technical support. Research was supported by the 
Army Research Office (W911NF-12-0607), the European 
Research Council and the Dutch Foundation for Funda¬ 
mental Research on Matter. E.K. was supported by a 
fellowship from the Nakajima Eoundation. 

I. Device details and additional data 

Device and measurement technique 

The device used for all the measurements shown in the 
main text, is based on a 12 nm wide Si quantum well 37 
nm below the surface in an undoped Si/SiGe heterostruc¬ 
ture, with two layers of electrostatic gates [see Eig.[5][a)]. 
Two accumulation gates (in yellow) are used to induce 
a two-dimensional electron gas ( 2 DEG) and a set of de¬ 
pletion gates (in gray), labeled 1-10 in Eig. [5][a), is used 
to form a single quantum dot in the 2 DEG (on the right 
side) with a charge sensor next to it (made by gates 4, 
5 and 6 ). Two 1 /im-wide, 200 nm-thick, and 1.5 fim- 
long Co magnets (in green), placed on top of accumula¬ 
tion gates (separated by an AI 2 O 3 layer), provide a stray 
magnetic field [see Eig. [5][b)]. The sample is thermally 
anchored to the mixing chamber of a dilution refrigerator 
with base temperature 25 mK and the electron tempera¬ 
ture estimated from transport measurements is 150 mK. 

Microwave excitation, applied to gate 6 in this ex¬ 
periment, generates an a.c. electric field, which 

makes the electron oscillate back and forth in the dot. 
Due to the gradient in the transverse magnetic field, 
dB_\_/dx^ estimated to be ^0.3 mT/nm [Eig. [5][b)], the 



FIG. 5. (a) False-color device image showing a fabricated pat¬ 
tern of split gates, labeled 1-10. For this experiment we create 
a single quantum dot (estimated location indicated by a blue 
circle) and a sensing dot (gates 4 and 5). The current through 
the charge sensor is recorded in real time for a fixed voltage 
bias of 500 fieV. The voltage pulses and microwave excitation 
are applied to gate 2 and 6 respectively. Green semitrans¬ 
parent rectangles show the position of two 200-nm-thick Co 
micromagnets. The yellow-shaded areas show the location of 
two accumulation gates, one for the reservoirs and the other 
for the double quantum dot region, (b) Schematic side view 
of the stray magnetic field generated by a single micromag¬ 
net, completely magnetized along the x-axis; B± represents 
the component of the stray magnetic field perpendicular to 
the static field. The external magnetic field Bext is applied 
along the x-axis [blue arrow in both (a) and (b) panels]. 

electron is then subject to an oscillating magnetic field 
^ \dB^/d.\ ^ eE^.^.\dB^/d.\ ^dicular 

^orb ^orb 

to the external static magnetic field. Here oc 
is the orbital level spacing, with lorb the typical QD di¬ 
mension. We notice that the amplitude of Ba.c. is propor¬ 
tional to the magnitude of the magnetic field gradient, to 
and to the amplitude of Ea.c.- 

We tune the right dot to the few-electron regime (with 
the left dot empty) and adjust the tunnel rate between 
the dot and the reservoir to be around kHz, so that 
dot-reservoir tunnel events can be monitored in real time 
by collecting the charge sensor current, read through a 
room temperature I-V converter, after a low-pass-filter 
with corner frequency of 30 kHz. The single-shot data 
[13,0 are processed on the fly by using an FPGA which 
directly counts the number of spin excited state events 
by using a threshold detection scheme. 

Estimation of the ratio of Rabi frequencies from CW 
measurements 

A typical CW spin resonance measurement is reported 
in Fig. El which shows the spin excited state probability 
as a function of the applied microwave frequency, /mw, 
and time. We note that the fluctuations of the two spin 
resonance frequencies, and are highly corre¬ 

lated; this suggest that the two resonances are affected 
by the same low-frequency source of noise (very likely 
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FIG. 6. The raw data on which Fig. 1(d) is based. Measured 
spin-up probability, as a function of applied microwave 
drive frequency fjviw and time (external field Bext= 560.783 
mT, power P = -12 dBm, microwave pulse duration 700 //s). 
Each horizontal scan in the figure takes ~2 minutes (200 cy¬ 
cles, in which each data point takes 2s), and the scan is re¬ 
peated ~60 times. The time reported on the y-axis is counted 
from the start of the measurement. 


hyperfine fluctuations) on the ^minute timescale. The 
trace for the second harmonic in Fig. 1(d) in the main 
text is obtained by averaging the data in Fig. [6] over the 
time axis. 

As reported in the supplementary information of 
Ref. [13 , the ratio of the steady-state spin-flip probability 
measured for the two valley states converges at low driv¬ 
ing power to r\{2ph)^ with rR(2ph) the ratio of the Rabi 
frequencies between two valleys [rR{2ph) = 

The ratio of the measured peak amplitudes in Fig. 1(d), 
which we write as A('L’+)/A(n_), is the product of 
r\(2ph) and the ratio of the respective populations at 
the end of the initialization stage, 


A{v-) 


9 

rR{2ph) X 


( 2 ) 


From a Gaussian fit to the spin resonance response for 
the two valleys in Fig. 1(d), we extract A{vj^)/A{ v-) ^ 
0.35. Furthermore, it is reasonable to assume that the 
ratio between the valley ground and excited state pop¬ 
ulations after the initialization stage, is the same when 
driving via the second harmonic as when driving via the 
fundamental. Therefore, we can use the ratio extracted 
for the fundamental in 1271 : 


g(4) _ g(2) ^ 0.3 ±0.1 
^ ^ ^ “ 0.7 ±0.1 
Prom eq. |2] and eq. [3] we obtain 


0.42 ±0.20. 


rR{2ph) = /f V/F = 0.90 ± 0.21. 


(3) 

(4) 


We remark that we also assume that the ratio of the 
Rabi frequencies between two valleys is the same at high 
microwave power as it is at low microwave power, which 
seems reasonable. 


Chevron pattern using second harmonic driving 


In Fig.[7l^a) we report the spin excited state probability 
oscillations as a function of the microwave burst time and 
detuning frequency (driving with Pmw= H-O dBm at the 
source). The quality of the data (stability of the measure¬ 
ment) is not high enough to extract independently the 
Rabi frequencies for the two valley states directly from 
the superimposed Chevron patterns (as was done in [27 1 
for driving via the fundamental harmonic). However, 
using the Rabi frequency ratio rR{2ph) ^ 1 extracted 
above, the ratio of initial populations of the two valleys 
of ^0.3/0.7 discussed above as well, and a Rabi 

frequency of 1.05 MHz for the valley ground state [esti¬ 
mated from a FFT along the MW burst time of Fig. [7l[a)] , 
we can simulate the two superimposed Chevron patterns, 
see Fig. HTb). This simulation agrees qualitatively with 
the data of Fig. (iKa), in particular it captures the slight 
asymmetry along the detuning axis, and the fact that 
mostly one Chevron pattern is visible. 


T2 estimation from Ramsey envelope decay and 
Ramsey simulation 

In order to get an alternative estimation of the dephas¬ 
ing time T|, we can perform a Ramsey measurement [see 
Fig. IHa) of the main text] and record the amplitude of 
the Ramsey oscillations as a function of the waiting time 
r between the two 7r/2 pulses, keeping /mw ~ /d^^- We 
show this data in Fig. [51 with the blue dotted curve rep¬ 
resenting the fitting relation = aexp[—(t/T 2 *)^] + c. 
From this fit we get a T 2 * of 780 ± 110 ns, compatible 
with what we estimated from the CW spin resonance 
linewidth in Fig. 1(d). On the same figure we report in 
black for comparison a fit with the exponential relation 
P^ = aexp(—t/T 2 *) + c. 


Echo decay of a qubit driven fundamental harmonic 
and comparison to results for driving at second 
harmonic 

In Fig. [9] we show a Hahn echo measurement real¬ 
ized driving spin resonance at the fundamental harmonic, 
recorded at the same magnetic field Bext as the measure¬ 
ment reported in Fig. 4(e) (driven by second harmonic). 
The T 2 ^^^{lph) extracted from a fit with a single expo¬ 
nential decay is similar to the T 2 ^^^{ 2 ph) extracted from 
Fig. 4(e). The decay obeys a single exponential similar 
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-1.5 0 -1.5 


detuning frequency (MHz) 


FIG. 7. (a) Measured spin-up probability, as a function of drive frequency Jmw and microwave burst time {Bext = 560.783 
mT, P = 11.0 dBm), (b) Simulated spin-up probability using population fractions 0.3:0.7, Rabi frequencies 1.0 MHz and 
1.1 MHz for the and v- valley states respectively, and the Larmor frequencies separated by ~1.04 MHz, extracted from 
low-power CW measurements in Fig. 1(d). 



t(hs) 


FIG. 8. Decay of Ramsey envelope as a function of the waiting 
time between the two Ramsey pulses. The blue dotted line is 
a fit of the data with the relation Pf = a exp[—(t/T 2 *)^] + c, 
with T 2 as a free parameter and c = 0.46, the average value 
for T T 2 (see Fig. |4|, which we also obtain keeping c a 
free parameter (c =0.46±0.09; T 2 *= 790±330 ns). The fact 
that the center of the oscillations is not 0.5 is attributed 
to initialization and measurement errors. 


to that observed in Ref. [l^, indicating the presence of 
a high-frequency decoherence process. 


n CO. 



t(hs) 

FIG. 9. Spin echo measurement for the ground valley-orbit 
state driven at the fundamental harmonic for Bext= 560.783 
mT and fMW = = 2/q^\ The observed times are 

similar for driving at the fundamental and second harmonic. 

II. Mechanisms mediating second harmonic 
generation. 

Here we elaborate on the three mechanisms discussed 
in the main text that can lead to second harmonic gen¬ 
eration. 

Position dependent magnetic field gradient 

From a simulation of the stray magnetic field of the 
micromagnet, see also the supplementary information of 
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Ref. [13, we can put an upper bound on the second 
derivative of the stray magnetic field with respect to the 
dot coordinates of cPB_\_/dx‘^ < 1 /iT/(nm)^. This is far 
too small to produce second harmonic Rabi frequencies 
that are comparable to the fundamental harmonic Rabi 
frequencies (which is what we observed experimentally). 

Tilted drive field 

If the a.c. magnetic field that is driving spin reso¬ 
nance is not strictly perpendicular to the static mag netic 
field, second harmonic driving becomes possible [33|. Eq. 
18 of Ref. [HI expresses the two-photon Rabi frequency, 
as a function of the drive strength uji oc Ba.c. 
and the drive frequency oo oc Btot^ for the case of an an¬ 
gle d between static and oscillating field of 45 degrees. 
They found ~ 2uj\Iuj. In our system the magnetic 

field gradient (which will provide Ba.c.{t) in combination 
with the microwave electric field) has components in all 
three directions of space. More specifically, from we 
have dBj j jdx ~ 0.2 mT/nm and dB_\_/dx ~ 0.3 mT/nm, 
from which we can estimate that 'd ^ 56 degrees. Fur¬ 
thermore, from we also get cji ~ 3 MHz and cj ~ 13 
GHz from which = 2ujiluj « 10“^, too 

small to explain the measured 

highest microwave driving field used in the experiment 
[see Fig. 2(b)]. 

Anharmonic confining potential 

Another possible mechanism which can give rice to the 
two-photon process in the Si/SiGe micromagnet-EDSR 
experiments can be related to the presence of an inter¬ 
mediate level (valley excited state spin down) lying be¬ 
tween the initial and final states (valley ground spin down 
and up) which can mediate the resonance process [pri¬ 
vate communication A. Palyi and M. Rudner]. If this 
is the case, by measuring the B-field dependence of the 
two-photon Rabi frequency, one should observe a non¬ 
monotonic f‘^Jli{Bext) dependence: the Rabi frequency 
should grow as half of the Larmor frequency approaches 
the valley splitting either from above or from below. 
Even if we did not perform a systematic study of the 
fiiabi ^ext 5 Rom the measurement we recorded we can 
exclude this striking dependence. 

In the absence of detailed knowledge of the shape of 
the confining potential, it is difficult to quantitatively 
estimate the magnitude of this effect. However, what 
we can say is that the details of the confining potential 
strongly influence the efficiency of second harmonic driv¬ 
ing: After tuning the device into a new gate voltage con¬ 
figuration (for which the valley splitting is higher than it 
is here), it has not been possible to observe again clear 


signatures of second harmonic driving. The new settings 
will likely also shift the average dot position in the mi¬ 
cro magnet stray field, thereby altering the tilt of the 
effective driving field as well as the position dependent 
field gradient. However, as discussed in the preceding 
sections, these two effects are several orders of magni¬ 
tude too small to explain the observed second harmonic 
driving. If the degree of anharmonicity of the confining 
potential dominates second harmonic driving, control of 
the anharmonicity will be required to make routine use 
of second harmonic driving in future spin qubit experi¬ 
ments. 
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